arXiv:1507.00752v3 [physics.optics] 16 Dec 2015 


Exact dipolar moments of a localized electric current distribution 


Ivan Fernandez-CorbatorQ 

Institute of Nanotechnology, Karlsruhe Institute of Technology, 76021 Karlsruhe, Germany 


Stefan Nanz and Rasoul Alaee 

Institut fiir Theoretische Festkorperphysik, Karlsruhe Institute of Technology, 76131 Karlsruhe, Germany 


Carsten Rockstuhl 

Institut fiir Theoretische Festkorperphysik, Karlsruhe Institute of Technology, 76131 Karlsruhe, Germany and 
Institute of Nanotechnology, Karlsruhe Institute of Technology, 76021 Karlsruhe, Germany 

The multipolar decomposition of current distributions is used in many branches of physics. Here, 
we obtain new exact expressions for the dipolar moments of a localized electric current distribution. 

The typical iirtegrals for the dipole moments of electromagnetically small sources are recovered as 
the lowest order terms of the new expressions in a series expansion with respect to the size of the 
source. All the higher order terms can be easily obtained. We also provide exact and approximated 
expressions for dipoles that radiate a definite polarization handedness (helicity). Formally, the new 
exact expressions are only marginally more complex than their lowest order approximations. 


The multipolar decomposition of a spatially confined 
electromagnetic source distribution is a basic tool in both 
classical and quantum electrodynamics [Il-@. On the one 
hand, the multipolar coefficients determine the coupling 
of the source to external electromagnetic fields. This is 
used in the study of molecular, atomic, and nuclear elec¬ 
tromagnetic interactions. On the other hand, there is a 
one-to-one correspondence between the multipolar com¬ 
ponents of the source and the multipolar fields radiated 
by it. This is exploited in the understanding and design 
of radiating systems. For example, in nanophotonics, the 
multipole moments of induced current distributions are 
used to study optical nano-antennas and meta-atoms [B- 
. The multipolar decomposition can be done in different 
ways, e.g. a Chap. 9] and [13, App. B, §4], resulting in 
integral expressions for the multipolar coefficients. The 
exact expressions are considerably simplified in the limit 
of electromagnetically small sources, but artificial scat- 
terers at optical frequencies are typically large enough to 
compromise the accuracy of the approximation. 


give the dipolar moments of electromagnetically small 
sources. The additional complexity is the appearance 
of spherical Bessel functions. We identify the spherical 
Bessel functions as the elements that perform the neces¬ 
sary selection of the appropriate momentum shell. When 
the spherical Bessel functions are expanded around zero, 
the typical approximations for the magnetic and electric 
moments of electromagnetically small sources are recov¬ 
ered as the lowest order terms in the expansion. The 
toroidal dipole is recovered as the second term in the 
electric case. All higher order corrections are easily ob¬ 
tained as successive terms of the expansions. We include 
integral expressions for the magnetic corrections of or¬ 
der and the electric/toroidal corrections of order k'^. 
We also provide exact and approximated expressions for 
dipoles that radiate a definite polarization handedness 
(helicity) [Eq. (HUl) and Eq. (HTl) ]. 

II. PROBLEM SETTING 


I. OUTLINE 

In this article, we obtain new exact expressions for 
the source dipolar moments [Eqs. (I20|) - (l22l) ]. In partic¬ 
ular, they are valid for any source size. We start our 
derivation in momentum space exploiting the fact that 
the fields radiated by the source at a given frequency uj 
are determined solely by its momentum components in a 
spherical shell of radius w/c, where c is the speed of light 
in the medium. We first obtain hybrid integrals in mo¬ 
mentum and coordinate space for all multipolar orders. 
In the dipolar case, we bring them to a form that is only 
marginally more complex than the typical integrals that 


We start by considering an electric current density dis¬ 
tribution J(r, t) embedded in an infinite, isotropic, and 
homogeneous medium characterized by real valued per¬ 
mittivity e and permeability p,. We assume J(r,t) to be 
confined in space so that J(r,t) = 0 for |r| > R. We 
consider its energy-momentum Fourier representation 


J(r,t) = TZ 

= n 


duj 

lo+ 


exp 

(Pp 


duj , ^ /■ ^ 

/ - 7 ^ exp i-iujt) / Jc.(p) exp (ip • r) 

Jo+ l-K J 


( 1 ) 
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and treat each uj term separately. The frequency ui and 
the three components of the momentum vector p are real 
numbers. The lower limit of the integral in duj excludes 
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the static case w = 0, which we do not treat in this paper. 
At each frequency w, the transverse electromagnetic fields 
outside the source are solely determined by the part of 
J^(p) in the domain that satisfies |pj = w/c. This result 
was obtained by Devaney and Wolf [^. We provide an 
alternative proof in App. 

We denote by Jt.j(p) the components of Ji.j(p) in the 
spherical shell of radius |p| = w/c. The symbol p rep¬ 
resents the angular part of the momentum vector p, i.e., 
the solid angle in the spherical shell. As usual, we define 
k = uj/c. 

We will expand Ji.c;(p) in an orthonormal basis for func¬ 
tions defined in a spherical shell: The three families of 
multipolar functions in momentum space [ 1 , B/.3] 


^jm (p) 


Zjm(p) 

W,™(p) 


1 

VJir+T) 


Ldjm (p) I 


ip X Xjm{p), 

pYjmip). 


( 2 ) 


The Yjm{p) are the spherical harmonics and the three 
components of the vector L are the angular momentum 
operators for scalar functions. 

Each of the vector multipolar functions in the three 
families is an eigenstate of the total angular momentum 
squared and the angular momentum along one axis q, 
for which we choose q = z. With Qjm(p) standing for 
any of the {Xj>„(p) 

; '^jm (p),Wj„(p)}: 

J^Qjm(p) =j{j + l)Qim(p), JzQjmip) = TOQjm(p), 

( 3 ) 

where j and m are integers, and m = For 

Xj„i(p) and Zjmip), j takes integer values in j > 0, 
while for Wj>„(p), j = 0 is also possible. 

The functions in Eq. are also eigenstates of the 
parity operatoi0: 


nx,„,(p) = -x,„,(-p) = (-i)^+ix,„,(p), 

nZjm(p) = —Zjmi — p) = (~l)'^Zjm(p), (4) 

nw,^(p) = -w,„,(-p) = (-i)^w,^(p). 

The polarization of Xjm(p) and Zjmip) is transverse 
(orthogonal) to p, and the polarization of Wj>„(p) is 
longitudinal (parallel) to p, as depicted in Fig. [T] In 
coordinate (r) space, this distinction corresponds to the 
distinction between divergence free (transverse) and curl 
free (longitudinal) fields. 

With the scalar product 

{A\B) = J dp Ai(p)B(p), (5) 


^ Their eigenvalues can be deduced from the parity transformation 
properties of a vector field in momentum space, i.e. nF(p) = 
—F(—p), and those of the spherical harmonics, ny^„i(p) = 
p) = (— 


where ^ denotes hermitian transpose, and p runs over the 
entire spherical shell, the three families together form an 
orthonormal basis for functions defined on any spherical 
shell in momentum space. 

We expand Jw(p) in this basis: 

j<v(p) = a^^J^jmip) + UHm^imip) + C^m^jniip), 

jm 

( 6 ) 

where, with standing for any of the 6“^, c“^}, 

<m = (QjmlL) = J dpQ]^(p)Jc.(p). (7) 


W(p) 



FIG. 1. The electromagnetic field radiated by a confined 
monochromatic current density Ju;(r) with Fourier trans¬ 
form Ji.,(p) only depends on the components of Ju;(p) 
in a spherical shell of radius |p| = oj/c. The relevant 
part of Ju;(p) can hence be expressed as a linear combi¬ 
nation of the momentum space vector multipolar functions 
{Xj7,.i(p), Zjm iP), (p)}, which form an orthonormal ba¬ 
sis for functions defined on the shell. The polarization vec¬ 
tors of Xjmip) and Zjmip) ^re tangential to the surface of 
the shell, i.e., orthogonal (transverse) to the momentum vec¬ 
tor p. The polarization vector of Wj„i(p) is normal to the 
surface of the shell, i.e., parallel (longitudinal) to p. 


The {a‘jmT^‘jm’djm} coefficients contain all the infor¬ 
mation about J(.j(p) so they must also contain all the 
information about the fields produced by it. As shown in 
[HI , the {ajm, b‘jm} determine the transverse electromag¬ 
netic field radiated by the sources at frequency w outside 
a spherical volume enclosing them: They are the coeffi¬ 
cients of the expansion of the transverse fields in outgo¬ 
ing electric and magnetic multipoles, respectively [2, Eq. 
9.122]. Therefore, the transverse components of Jt.j(p) 
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determine the transverse components of the electromag¬ 
netic field at frequency w outside the source region. The 
longitudinal electric field with |p| = w/c is zero outside 
the source region. While the longitudinal degrees of free¬ 
dom of i.e. the are not necessarily equal 

to zero, the field that they generate outside the source 
region is canceled by the fiel d g enerated by the charge 
density. This can be seen in |ld §13.3 pl875-1877], and 
in [ 13 , App. C] where the cancellation is shown to be a 
consequence of the continuity equation. We will keep the 
in the discussion both for completeness and because 
they play an important role in understanding the split of 
the into electrical and toroidal parts |14l-[lq. which 
we discuss in [II- 

The {a“^, b'j^} coefficients are a valuable source of in¬ 
formation in many branches of physics. In molecular, 
atomic and nuclear physics, the coefficients 

are used to describe the interaction of systems of charges 
with external electromagnetic fields, e.g. [1, Chap. 10], 
d, IV.C.2c)] and [1, Chap. 7]. In classical electrodynam¬ 
ics they are used to describe radiation by source distri¬ 
butions, e.g. [I|, Chap. 9] and [1, Chap. 9]. In nanopho¬ 
tonics, they are used to study and design the response of 
individual artificial nanostructures. 

Given Jij(r), there exist exact expressions for the 
coordinate space integrals, e.g. [1, Eq. 

(7.20)] 


Ajm = -^ / (V X jj{kr)Xj^{r))^ 3^{r), 
= J d^r jj (fcr)X]^(f)J<^(r), 


( 8 ) 


^10 

u; 

^1 — 1 


(fir Ja;(r), 


^10 

^1-1 


d^r r X J<.,(r), 


( 10 ) 

where we have chosen the spherical vector basis. We will 
work in this basis throughout the article. Appendix [C] 
contains auxiliary expressions. 


III. EXACT DIPOLAR MOMENTS 

We will now obtain exact expressions for the 
dipolar vectors [a^i, o^q, [^n, & 10 ) ^“- 1 ]^ and 

[cii, c^Q, as coordinate space integrals of functions 

of Jcj(r). While these expressions are, as Eq. (|S]) and Eq. 
(P|). valid for any source size, they are only marginally 
more complex than their kR <C 1 limits: Namely, they 
contain spherical Bessel functions. As far as we know, 
these expressions have not been reported before. 

We start from Eq. 0, where we substitute 

to get 

Q]m(p) J d^r J^(r) exp (-*^P • r) ■ 

( 12 ) 

The condition jpj = w/c is enforced in the argument of 
the exponential. We now substitute the exponential for 
its expansion in spherical harmonics 


or, a Eq. (9.165) without the magnetization current 
therein], 

Sjm = / d^r jj{kr)Y*^{r)L ■ (V x J<^(r)), 

bjm = /.fi ... / d^r jj{kr)Y*^{r)L ■ J<^(r), 

(9) 

where the tildes and carets in the left hand sides indicate 
different normalizations, r = jrj, f = r/jrj is the angular 
part of r, and k = uj/c throughout the article. 

The expressions in Eq. ([5]) and Eq. m are valid for any 
source radius R. Eor electromagnetically small sources 
where kR ^ I, they can be reduced to the simpler well 
known expressions that are obtained in a Chap. 9] and 
a Chap. 9] by starting with the equation for the vector 
potential as a function of (r) in the Lorentz gauge [Eq. 

(IA21) ]. and expanding in powers of fcjr — r'j. 

Eor example, when kR <C 1 the source dependent terms 
of the electric and magnetic dipole moments are 


exp (-*^P • r) = (47r)^(-i)'y4(f)y|^(p)j^(fclrl), 
Z,m 

(13) 

where jY(-) is the Tth order spherical Bessel function of 
the first kind. The result is: 

vW . ^ 

J dp Q]m(p)lrm(p) J d^r J„(r)r4(f)j;-(A:r). 

Ini 

(14) 

Equation (HI is an exact expression for the 
djm} Coefficients in terms of integrals in both 
momentum (shaded area) and coordinate space. As 
shown in App. [B] only terms with / = j contribute to the 
6“^, while the and get contributions from both 
I = j — 1 and I = j -I- I. Additionally, it is possible to fur¬ 
ther simplify Eq. (imi in the dipolar (j = I) case without 
making any approximation. We now present the deriva¬ 
tions for the magnetic dipole Appendix |D| contains 
the derivations for and H also contains the Cqo 
case. 
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For the magnetic dipole, we particularize Eq. (1141) for 
Qjmip) Xjm(p) and j = 1, which implies 1 = 1: 


. \/(27r)= 
* 47r 

Tn—1 


— 


dp X|„(p)Fi7n:(p) J cfr 3 ^{r)Y*^{r)ji{kr). 

(15) 

Explicit expressions of Xim(p) can be obtained using 
Eq. (IC6I) and then used to write the momentum integrals 
in the shaded area of Eq. (13 as 

n t 


m = 1 —>■ 


m = 0 


m = —1 — 


dp 


yio(p) ' 

V2 

yidp) 

V2 

0 


hlm,(p), 


/ dp 

1- 

1 >o 

1 

1_ 

J 

yll(p) 

L V2 J 


n t 


Elm(p), 


(16) 



r 0 1 

y dp 

Ti-i(p) 

V 2 

yio(p) 

L ^ J 


n t 


yirniP), 


which can be easily solved for each m € {—1,0,1} using 
the orthonormality properties of the spherical harmonics: 
/ dp Yimip)Yj^{p) = Srnmh- They result in three vec- 
tors for each m case, which we list here as row vectors. 
From top to bottom, the three row vectors correspond to 
m = 1, 0, —1: 


Considering the expression for the cross product in spher¬ 
ical coordinates [Eq. (|C4p ]. we can finally write Eq. (fTg|) 
as: 


uii 

^10 

Lu; 

Ol_l 


= -^J d^r f X J^(r)ji(fcr). (20) 


The expressions for [«ii: “ioj 
[ cii, c^Q, can be obtained by similar, although 

more involved, procedures. We provide the derivations 
in App. [D] The results read: 


Oil 

Uio 

n ^ 
“l-l 


1 


rV3 . 


d^r Ja;(r)jo(A:r) 


z=o 


2^-^/ d?Y |3 Cj„(r) f-J<^(r)| j2(fcr) 


1^2 


( 21 ) 


and 


-10 


r-\/6 . 


d^r J^(r)jo(fcr) 


!=0 


|3 f-J,^(r)|}2(fc?’) 


1=2 


( 22 ) 


TO = 1 —>■ 


TO = 0 —>■ 


TO = —1 —> 


1 

71 

1 

71 

1 

71 


( 

0 

1 

0 

) 

( 

-1 

0 

0 

), 

( 

0 

0 

0 

) 

( 

0 

0 

1 

) 

( 

0 

0 

0 

), 

( 

-1 

0 

0 

) 

( 

0 

0 

0 

) 

( 

0 

0 

1 

)• 

( 

0 - 

-1 

0 

) 


(17) 


Having solved the momentum space integrals in the 
shaded area of Eq. (USD, the summation in m can now be 
done. With Jcj(r) = [J^, Jq , ^ and, as in Eq. (|C5p . 


f = 


r 



Y,\{v) 

n*o(f) 

Ti*-i(f) 


(18) 


the result of the sum reads 

^11 = d^Y {J^f^ - J^ro) ji(fcr), 

&10 = ^ y d^Y ji{kr), ( 19 ) 

^1-1 = {J-iro - Jof-i)ji{kr). 


where the contributions coming from / = } — 1 = 0 and 
{ = j + 1 = 2 are indicated. The dot product fUtj(r) is 
simply equal to f^J„(r) in Cartesian coordinate^ 
Equation (1^ . Eq. (1^ . and Eq. (1^ are exact. In 
particular they apply to a source distribution of any size. 
They are also simpler than the corresponding exact ex¬ 
pressions obtained from Eq. ([3 or Eq. We note 

that Eqs. (1201) to Eq. (1^^ should also be reachable from 
the coordinate space integrals of Eq. (jS]) or Eq. ® . Our 
route through momentum space explicitly exploits that 
the contributions to the only come from the Eourier 
components of the source in the domain |p| = w/c. This 
restriction is imposed in the exponential of Eq. m and 
determines the argument of the spherical Bessel func¬ 
tions ji{kr) in Eq. (dH), which then appear in Eqs. ([^U)) . 


^ This can be seen using the complex unitary matrix that 
transforms spherical into Cartesian components, written on 
the right of Eq. Elli. Calling such matrix C, we start 
with the dot product in spherical coordinates and transform 
it to Cartesian coordinates using that (C^C) is the identity: 
ftJij(r) = ft(C^‘C)Ju;(r) = (Cf)t(CJuj(r)) = rl^^^Jui(r)cart = 
^'cart^(^ (r)car£, where the last equality follows because the Carte¬ 
sian coordinates {x,y,z) are real. 
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Eq. (1^ . and Eq. (|2^ . We can deduce that the spheri¬ 
cal Bessel functions must be responsible for rejecting the 
IpI ^ w/c components present in Ji^(r). We now provide 
a more formal proof of their role. 

In the expression of in Eq. m, the dependence 
on the current density is contained in the integrals 


(fr 3^{r)Yi*^{r)ji(kr). 


(23) 


the same power of k we obtain: 

■ 

Oil 
^10 
Oi_i 




Z=0 


xii 

-.UJ 

ho 


^1-1 


2ttV3 


k J (Py r X Jaj(r), 


(27) 

(28) 


We then write Jaj(r) as an inverse Fourier transform and 
expand its exponential exp (ip • r) as in Eq. (IISl), except 
that now |p| is not restricted to Lojc. After rearranging 
the integrals we get: 



(24) 

The shaded d^r integral can be solved by splitting it 
into its radial and angular parts (/ d?Y = dr f dr'). 
First, the angular part is solved through the orthonor¬ 
mality of the spherical harmonics, which forces {l,rn) = 
(l,m). The remaining radial integral has a formal solu¬ 
tion as a radial Dirac delta distribution [l3, Eq. (4.1)] 

j drr‘^ji{\p\r)ji{kr) = ^S{\p\ - k), (25) 

which enforces the jpj = k = ojjc restriction in Eq. (l24l) . 
namely: 

--^=J d^pJ^(p)Y;^(p)^S{\p\-k) = 

y dpYil.ip) dp p^3^{p)^S{\p\-k) = 


I r-2r2j,^(r)}, (29) 

r^oj^2 

J^(r) (30) 

X ^ > 

r-0 

^ ^ ■ (31) 

V 

Equation (07)) . Eq. (05)) and Eq. (051) are, respec¬ 
tively, the well known approximated magnetic, electric, 
and toroidal dipole moments of electromagnetically small 
current distributions. We note that the electric dipole 
contains contributions only from I = 0 while the toroidal 
dipole has contributions from I = 0 and 1 = 2. 

The small argument approximation causes two kinds of 
inaccuracies. On the one hand, entire integral terms are 
neglected. For example, the toroidal term in Eq. (155) dis¬ 
appears in a lowest order approximation. On the other 
hand, some components with jpj ^ ujjc will leak into 
the dipole moments. This happens because the approxi¬ 
mated expressions of the spherical Bessel functions do not 
correspond to momentum space Dirac deltas d(|p| — k). 

Approximations with increasing accuracy are obtained 
in a straightforward way from the exact Eqs. (1551) to Eq. 
(1551) . It is a matter of taking more terms in the expan¬ 
sions of the spherical Bessel functions. For example, the 
(/cr)^ correction to Eq. (1571) reads 


'-10 

" 1-1 


(26) 


The ji{kr) functions from Eq. (1141) find their way into 
Eq. (1551) , and become one of the pieces needed to obtain 
the Dirac delta (5(|p| — k) which filters out the jpj 7 ^ uj/c 
components of Ja;(i’). 



the {kr)^ correction to the total in Eqs. (158l) - (1551) 
reads 

(33) 


IV. ELECTROMAGNETICALLY SMALL 
SOURCE APPROXIMATION WITH 
INCREASING AGGURACY 

We now make the small argument approximation to 
the spherical Bessel functions in Eqs. (isni-diii) and keep 
terms up to second order: jo{kr) Ri 1 — (fcr)^/ 6 , ji(kr) k. 
kr/Z and j 2 (kr) k. (fcr)^/I5. After grouping terms with 


and the [kr)'^ correction to the total in Eqs. (1551) - 
(1511) reads 

|[rtj^(r)]r+ij,.(r)r2|r2. (34) 

The above corrections to and coincide up to 
normalization factors with the mean square radii in [isl . 
App. C], where they are derived in a different way. 
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We now use our results to compute the magnetic 
dipole moment of a current distribution with a previ¬ 
ously known analytical solution, verify that the result 
coincides, and compare it with two approximated solu¬ 
tions for electromagnetically small sources obtained from 
taking the first and the two first terms in the expansion 
of the spherical Bessel functions. 


V. EXAMPLE 


Let us consider an infinitesimally thin circular loop of 
current with implicit time dependence exp(—iwt). The 
loop has radius a and lies on the plane perpendicular to 
the z axis (see the inset in Fig. The expression for 
its current in spherical coordinates is 

3u,{r) = ^IoS{r - a)-6{6 - (35) 

r Z 

where cj) = [—sin </>, cos (^, 0]^, (j) = arctan(^) and 6 = 
arccos(^). 

The exact value of its magnetic dipole moment is ob¬ 
tained after calculating the integral in Eq. (EUl) : 

m = 2,V3Ioaji{ka). (36) 


We obtain a first small source approximation by us¬ 
ing Eq. (07)) and a more accurate second one using the 
incremental correction in Eq. OSH) 


mfea«i = zv^7o —, 

mia«i = [l - [kaf/lQ] . 


(37) 


The same results are obtained by taking terms up to ka 
and (fca)^, respectively, in the Taylor series of ji(fca) in 
Eq. (1551) . This latter approach relies on the existence of 
an exact closed form solution and is hence not general. 

The exact value of Eq. (1551) coincides with the one 
calculated in [I^, §13.3 pl881] up to a numerical factor 
that can be traced back to a different normalization. In 
this simple example, the relative error incurred due to 
the small source approximations is equal to the relative 
error incurred when approximating the first order spher¬ 
ical Bessel function. Figure [5] shows the relative errors 
incurred when taking only the first term in the expan¬ 
sion [ji(fea) ~ ka/S\ and when taking the first two terms 
{ji{ka) ~ (ka/i) x [l — (/ca)^/10] }. We see that, if we 
take only one term, a 10% relative error is incurred when 


the diameter of the loop is approximately 30% of the 
wavelength. When taking two terms, the 10% relative 
error is reached when the diameter is approximately 70% 
of the wavelength. We note that in this example the cur¬ 
rent is concentrated in the most exterior region of the 
object. When this is not the case, e.g. in a homogeneous 
current distribution within a sphere of diameter 2a, the 
relative errors should be smaller. 



FIG. 2. Relative error in the magnetic dipole moment of an 
infinitesimally thin circular current loop of radius a (shown 
in the inset) due to the small 27ra/Ao approximation. Solid 
red line: Error due to taking only the first term in the small 
argument expansion of the spherical Bessel function in Eq. 
dsni). Such first order gives the typical integral for the mag¬ 
netic dipole moment of electromagnetically small sources [see 
Eq. 0271)]. Dashed black line: Error due to taking the first 
two terms in the expansion, i.e. Eq. (1271) plus Eq. (13211 . 


VI. RESULTS FOR HELICITY MULTIPOLES 

There is some recent interest in the use of helicity for 
the study of interactions between matter and electromag¬ 
netic fields [Tsl - f^ . Due to its fundamental relationship 
with electromagnetic duality, the helicity formalism is 
also very useful when discussing dual symmetric systems 
[lEllJ, e.g. Huygens surfaces [H, |2^. We now extend 
our results to the dipoles of well defined helicity. 

Multipoles of well defined helicity are an alternative to 
the multipoles of well defined parity. The two sets are 
related by a change of basis, which we write for both the 
coefficients and the Qjm(p) functions: 
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UJ jm ' ^jm ^ 

r'+ 


yjm+ ^ 


V2 

to jrn ^ 

G-^(P) 

^jm(p) ^jm(p) 

9jm- ^ 


U) _ OJ , 

yjmO ^j'm 

G°„(p) 

— 'Wjj7j(p) — 


1 -|- ipx TjY^jui 

\/j(j + 1)’ 

1 ipx 

V2 \/j(/+1)’ 


(38) 


where 1 is the 3x3 unit matrix. 

The G^^(p) in Eq. (1551) and the Qjm(p) have the 
same properties under rotations. They differ in their par¬ 
ity and polarization properties. Instead of eigenstates of 
parity, the G^^(p) are eigenstates of the helicity opera¬ 
tor with eigenvalue A. This is obvious from the rightmost 
expressions in Eq. (1551) since the helicity operator A in 
the momentum representation is ipx : 

A=^^->ipx, (39) 

where J and P are the angular and linear momentum 
vector operators, respectively. 

The two transverse families of this alternative basis, 
G^^(p), correspond to multipolar components g‘jm± that 
radiate fields of definite polarization handedness (helic¬ 
ity) A = ±1 [13, App. A]. 

The extension of our dipolar results to the helicity ba¬ 
sis is straightforward. According to the third line of Eq. 
(1551) . the result for A = 0 is Eq. (1551) . The exact expres¬ 
sions for the transverse dipoles with helicity A = ±1 can 
be obtained using Eq. (150l) . Eq. (15T1) and Eq. (1551) : 


The approximated expressions up to order are: 



VII. CONCLUSION AND FUTURE WORK 

In conclusion, we have obtained new exact expressions 
for the dipolar moments of a localized source distribu¬ 
tion. These expressions are simpler than the ones re¬ 
ported to date. They are only marginally more complex 
than the typical integrals for the dipole moments of elec- 
troinagnetically small sources and allow to easily obtain 
approximate expressions with increasing accuracy. Our 
results can be applied in the many areas where the dipole 
moments of electrical current sources are used. 

In future work, we aim to obtain new exact expressions 
for general j-polar order and use them in applications like 
for instance in the study of the scattering properties of 
nanostructures. 
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Appendix A: Fields produced by time varying 
sources: Only Fourier components with |p| = uj/c 
contribute 
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We consider a electric charge and current density dis¬ 
tributions p(r, t) and J(r, t) embedded in an isotropic and 
homogeneous medium with constant and real permittiv¬ 
ity e and permeability We assume them to be confined 
in space so that p(r, t) = 0 and J(r,t) = 0 for |r| > R. 
We consider the following Fourier decomposition: 


p(r,t) = 


= 7^ 


r 

0+ 

OO 7 

duj 


exp(-zu;t) p^(r) 

(Pp 


f duj , . f Pp 


J(r, t) = TZ 


r 

/o+ 


= 7^ 


exp(-iwt) J^(r) 


duj , , f d-^p ^ , 

/ - 7 = exp {-lujt) / exp (zp ■ r) 

Jo+ J ^{2'kR 


(Al) 


The lower limit of the integral in doj excludes the static 
case a; = 0 . 

Devaney and Wolf 0 proved that, outside the source 
region, the transverse parts of the electromagnetic field 
produced by the source at frequency uj are determined by 
the transverse components of Jc^(p) that meet |p| = w/c, 
where c = Ij is the speed of light in the medium. 

We now provide a different proof which uses the po¬ 
tentials instead of the fields and shows the selection of 
the IpI = uj/c components through the appearance of a 
radial delta distribution. We prove that the only parts 
of the sources that contribute to the scalar and vector 
potentials in the Lorenz gauge are those in the domain 
IpI = w/c. The electric and magnetic fields obtained from 
the potentials are hence also determined by the compo¬ 
nents in the momentum space shell with radius |p| = w/c, 
which means that the result is independent of the choice 
of gauge. 

In the Lorenz gauge, and with implicit monochromatic 
exp {—i(jjt) dependence, the sources in Eq. dSH) generate 
the following scalar and vector potentials: 


^ J d^r'p^{r') 


, exp {ik\r — r'|) 


47r|r — r'l 


A^{r) =p d^v' J^(r') 


, exp (ifc|r — r'l) 


(A2) 


47r|r — r'l 


where k = oj jc. 

Following Jackson’s steps, we use the expansion of 
exp {ik\r — r'|) /(47r|r — r'|) in [2, Eq. 9.98] 


= £ yut)YUt') 


47r|r — r' 


l=Q 


(A3) 


to get to a Eq. 9.11]: 

Ac„(r) = 

ipkY^hf\kr)Yim{v) j d^r' 3^{r')ji{kr')Yi*^{r'), 

l.m V ^ 

Tim 

(A4) 


where I and m are integers, and j/(-) are the Lth 

order spherical Hankel and Bessel functions, respectively, 
f = r/|r|, f' = r'/|r'|, r = |r|, r' = |r'|, and Yjg are the 
scalar spherical harmonics. 

Let us now consider the integral labeled as Tim in Eq. 
(IA4I) for a given term {l,m). We use the inverse Fourier 
transform of Ja;(r') 

exp(ip-rO, (A5) 

J V(27r)'^ 


and the expansion of the exponential exp (ip • r') in 
spherical harmonics, 


3 — 1 


exp (ip • r') = (47r) ^ ^ i^Yj^{r)Yi;^{p)ji{\p\\r'\), 

(A 6 ) 


m——l 


to get 

FiTn _ 
47r 


Trn V V / 

(A7) 


We stress that Jcj(p) in Eq. (IA5I) . and hence Ji.j(r') 
in Eq. (IA4II and Tim in Eq. (IA7D . may contain con¬ 
tributions from momenta p such that |p| ^ ijjjc. The 
following steps show that these contributions are filtered 
out and that At,j(r) depends only on the components of 
Ji.j(p) with IpI = w/c. 

We take Eq. (IA7I) . split the integral in PY into ra¬ 
dial and angular parts (/ PY = dY Y^ J df'), and 
solve the angular part through the orthonormality of the 
spherical harmonics / drYim{r)Yi^{r) = SuSmm- After 
this, the only term in the sum on I and m that does not 
vanish is the one meeting I = I and Wi = m\ 



(A 8 ) 


The crucial step is that the integral in the shaded box 
of Eq. (IA 8 I) has a formal solution as a radial Dirac delta 
distribution [TtI . Eq. 4.1]: 

J dY{Yyji{\p\Y)ji{kY) = ^S{\p\- k). (A9) 

This S{k — IpI) term discards all momenta contribu¬ 
tions from outside the spherical shell |p| = k = uj/c 
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in Eq. (IA 8 I) . To show it explicitly, we split the in¬ 
tegral in into radial {p = |p|) and angular parts 

if d^P = fr dpp"^! dp): 

^im = J dp Yi*^ip) j dpp^3^pp)^5(j)-k) = 
J dp J<^(p, IpI = fc)r,;,(p). 

(AlO) 

Since this conclusion holds for all values of {I, m) in Eq. 
(E4l) . it follows that the vector potential is completely 
determined by Jaj(p, |p| = k), i.e., the components of 
Jw(p) on the momentum shell of radius |p| = uj/c. 

The same conclusion is valid for the scalar potential 
in Eq. (IA2I) . This can be seen noting that none of 
the steps in the previous derivation needs the fact that 
Jtj(r) is a vector. The same steps can be taken for the 
scalar charge density p^ (r) which generates the scalar po¬ 
tential in Eq. (IA2I) . Regarding its inverse Fourier trans¬ 
form 




d^p 

'/Wf 


p<^(p)exp (ip - r), 


(All) 


the conclusion in this case is that (pui (r) only depends on 
the momentum components of the charge density /9i.j(p) 
in the momentum shell of radius |p| = lo/c. 

Since both scalar and vector potentials (pij(r), Atj(r)) 
depend only on the source Fourier components in the 
domain |p| = w/c, the same will be true for the electric 
and magnetic fields computed from them: 


E^(r) = ia;A,^(r)-V())c,;(r), B^(r) = VxA„(r). (A12) 


It is hence clear that the conclusion is gauge indepen¬ 
dent. It is also clear that the derivation applies to both 
transverse and longitudinal components of the electro¬ 
magnetic field, but the longitudinal electric field with 
IpI = uj/c is zero outside the source region. This can 
be seen in [I^, §13.3 pl875-1877], and in [I^ App. C], 
where the cancellation is shown to be due to the conti¬ 
nuity equation. 


Appendix B: One term in fcJA) two in and 

We show that, for and only terms with 

I = j — 1 or I = j + \ can be different from zero in 
Eq. (fT4l) . and that for 6 “^, only I = j contributes. For 
this, we will write the momentum space integrals in the 
shaded area of Eq. (fHl) as integrals of triple products of 
spherical harmonics. These integrals have an exact ex¬ 
pression involving a product of two 3j-Wigner symbols. 
The requirement that one of the 3j-Wigner symbols be 
non-null results in the aforementioned relationships be¬ 
tween j and 1. 


Particularizing the shaded integrals in Eq. (|T4)) to 
Wj>„(p) and ZjmiP) 

J dp Wj^(p)ty,^(p) = y dp [pYj^ip)]^ 

j dp Zl^{p)Yi„{p) = J dp [ip X Xj>„(p)]^ ^m(p), 

(Bl) 

we find that each of their three components contains 
either one [in the Wjm(p) case] or a sum of two 
[in the Zjm{P) case] triple products of spherical har¬ 
monics like Y*pY*^Yij., which can be also writteiJl as 
(^—I'^P+iYi-pYj-qYij,. The result of the integral of the 
product of three spherical harmonics is [H, p. 700] 

J dfl Y/jmi (f^)Ei2m2 (f^) = 

j {2li + l)(2l2 + l)(2i3 + 1) I'll I 2 /sN / li I 2 /a \ 

V 47r 1^0 0 Oy m 2 m 3 J ' 

(B2) 

where ( ) is the 3i-Wigner symbol. 

\qi 92 93/ J B J 

In our case, we see from the right hand sides of Eq. 

m that ji = 1 from p [see Eq. deg], j 2 = j from 

ZjmiP) or Wj>„(p), and js = I from Yf:^{p). We now 

consider some of the conditions for the first 3j-symbol in 

Eq. (|12]) 

G i 0 (‘=3) 

to be different than zero. Namely [ 2 ^ p. 1056]: 

lii - J 2 I < 33 < ji + h |1 - j| < i~< j + 1 (B4) 

which, if J > 0 restricts I to be j — I, j or j -|-1 and, when 
j = 0 in the longitudinal case, forces 1 = 1. Furthermore, 
because of the zeros in Eq. (IB3D , I -I- j -I- T must be an 
integer multiple of 2, which then forbids I = j when j > 0. 
All together we obtain for the 'Wjm (p) and Zjm (p) cases 
the restrictions: 


/ = j - I or J + 1 if J > 0, 
1 = 1 if j = 0 . 


(B5) 


In the Xj>„(p) case, the integrals in the three compo¬ 
nents of 

J dp X]^(p)r|^(p) (B 6 ) 


contain a product of two spherical harmonics, but the 
third one can always be assumed to be the constant I = 
•\/47rFoo- In this case the restriction of Eq. (IB4I) forces 
i=3- 

\0-3\<I<3 + 0 ^ I=j. (B7) 


3 Using that 










11 


Appendix C: Auxiliary expressions in the spherical 
vector basis 


in Eq. ([2]) and the expression of the angular mo¬ 
mentum vector operator L in spherical coordinates 


We write a vector a in the spherical vector basis as: 

a = ai©! -I- cto®o O— 1 ®—ij (d) 


with 


— Lx-\-iLy 


-^'down 

V2 


V2 

L, 

= 

-^0 

Lx~\-‘iLy 


-hup 

L \/2 J 


L V2 -1 


(C7) 


x-hiy 

©0 = 2 (C2) 

X- iy 

V2 ■ 

This choice of basis induces the following relationships 
between the Cartesian and spherical coordinates of a in 
the spherical and Cartesian basis: 


where Lup = + iLy and Tdown = — iLy are the 

angular momentum ladder operators 


f ^j{j -I- 1) - m(m -I- l)yj(^+i) if|m-|-l|<j 

\ 0 else 


-^'down^^'m — 


f VjU + 1) - - l)7S'(m-l) if |UI - 1| < j 

\ 0 else 

(C8) 


ai 


- -1 
V2 

V2 

O' 


ax 


a 

flo 

= 

0 

0 

1 


Qy 

5 

a 

0-1 


1 

LV2 

V2 

0 


_az_ 


a 


o 

_1 


Oi 

— n ^ 

V2 ^ V2 


Oo 

O 

1—1 

o 


_a-i_ 


(C3) 

In the spherical basis, the components of the cross prod¬ 
uct of two vectors are[3: 


a X b = i 


aibo — agbi 
aib-i — a-ibi 
O'ob-i — ci-ibo 


(C4) 


Let us now write some explicit expressions for p and 
Xjm(p) that we use in the text. 


Pi 


-Vi-1 

Fk 

' V* 
r 11 

Po 

= ^3 

no 

= 2V3 


P-i_ 

V ^ 

-Vii 

V 

J- 1_1_ 


1-1 

(C5) 


j(j + l)-m(m-l) /-^ 

2 ^j(m-l)(P) 


Xjm(p) 


1 

VW+T) 


mVjm(p) 


(C6) 

Equation (IC5I) follows Eq. (IC3I) . the expressions of 
Yim in Cartesian coordinates and the property Y*^ = 
{—lyYi-q. Equation (|C6p follows from the definition of 


Equation (El is obtained from the relations ei x e_i = zeo, 
ei X Go = zei and e_i x gq = —zg_i, which follow from Eq. 
del and the cross products between the Cartesian basis vectors 

{x,y,z}. 


Appendix D: Expression of selected tensors as 
spatial integrals 

1. Case CLXrn 

As shown in App. |B] only I = 0 and I = 2 can have 
non zero contributions to aim- That is 

a-im = a-J=^ + a-J=^. (Dl) 

We start with 1 = 0. From Eq. (HI, and since Yqo = 

\/ 


J=o 


1 


dpZ\^{p) (fr 3^(r)joikr). 

(D2) 


The explicit expressions for Zi^ = ip x Xi^ are 

27r 


m = l-)- ipx Xii(p) = — 


m = 0^ ipx Xio(p) = -\l Y 


Yit - >^1^1-1 


-YiiYio 

y2 

''^11 

■ YioYi_i ■ 
-2YiiYi_i 
YioYii 



[ Y2 1 

-^1-1 

Vt 

1 

o 

1 

nt - niVi-i 
(E 


The relationship YJ* = (—l)^y/_g, and the orthonormal¬ 
ity of the spherical harmonics allow us to solve the mo¬ 
mentum space integrals in the shaded area of Eq. (imi) . 
and immediately reach 

I=Q 

=-;^y J<^(r)jo(A:r). (D4) 


u^io 

UJ 

^1—1 
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In the I = 2 case 


VWf 1=2 


47r 

rn—2 


- E 


dV Z\^{p)Y2rn 


d^r 3^{r)Y*^j2ikr), 


(D5) 


the shaded momentum space integrals contain triple 
products of spherical harmonics and can be solved us¬ 
ing 


dpYi^rmip)Yi^ 

m2 iP)Yl3 m 3 (p) 

(2/1 -|- 1)(2Z2 -b 1)(2^3 -|- 1 ) //j I2 I3 

47r I 0 0 0 


h h h 

mi m 2 m 3 . 

(D 6 ) 


where 


is the 3j-Wigner symbol. 


Jl J2 J3 
92 932 

They result in five vectors for each m case, which we 
list here as row vectors. From top to bottom, the row 
vectors correspond to m = 2 , 1 , 0 , — 1 , — 2 : 

I 


1=2 


1 


47r 


m = 1 —>■ 


m = 0 


m = —1 —>■ 


(0 0 V6 ) 
( 0 -73 0 ) 
(10 0 ) 
( 0 0 0 ) 
( 0 0 0 ) 

( 0 0 0 ) 
(0 0 73 ) 
( 0 - 20 ) 
( 73 0 0 ) 
( 0 0 0 ) 

( 0 0 0 ) 
0 0 0 ) 
0 0 1 ) 


7^ 


-1 

7^ 


-1 


(D7) 


( 


( 


( 0 -73 0 ) 
( 76 0 0 ) 


The summation in m in Eq. (Ids!) can now be done. 
With J^(r) = [J-, Jp-, and Y^*^ = (-1^.^, it 

reads 


1^2 


^10 


1^2 


7(^7M. 

1 47r 

7(^7M. 

1 47r 

7(^7M. 




(76J^ir2_2 + V3J^Y2-i + J^Y2o) j2{kr), 
d^r (-73J“iT2-i - 2J^Y2o - V3J‘^Y2i) j2{kr), 
d^r (^J‘^iY2o + V3Ji^Y2i + 76JrT22) j2ikr). 

I- 


(D 8 ) 


We now use the following relationships: 
^22 = 


IOtt 


Yi\, Y 21 = 


207r 


FioTii 


r 2 o = 7^(y7-^ 


y _ 7 IOtt 2 y - y 

^2-2 — Y ^2-1 — Y -^riori-1, 

which we substitute in Eq. Idhi and get 


(D9) 


uj 1—2 
3-11 


73 . 


7r 


Yi-i (+ JqYiq^ + 


JUJ 

Jl 


— 
-^10 


Att 


1^2 


^10 


^ 1—1 


^ j A no + J,"y„) + .js (ri, - T) 




/ 1 


Yii {J^Yio + J^Yii) +^{Y,i- — 


hikr), 
hikr), 

j 2 {kr). 


(DIO) 
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The expressions in the shaded areas of Eq. (IDIOI) can 
be completed to YnJ'^ + Tio+ Yi-iusing terms 


to their right. In the case of the ^ the completion is iJ-coo = I (fr 2 ,/- + J^Y^^ - JiY^*_^) ji{kr) 


straightforward. For the other two cases one uses that 


47r 


= |yio|2+|Fii|2+|ii-il 


1 


1 


= / d^Y 


= :^+2riiri_i. 


f^Ja;(r) ji(fcr). 


47r 27r 


(Dll) 


(D17) 


Finally, noting that 


Fnjr + FioJo“ + Fi_i^V ^ , (012) 


The first term in a small kr expansion of cqo will be of 
order fc: 


Coo ~ —* 


we reach the final result 
n 1=2 


2 k 
TT 3 


J dh [r^J^(r)] . (D18) 


a 


a 


1-1 


1 


271^3 , 


J d^r |3 f'''Ja;(r) r - j 2 {kr). 


(D13) 

The sum of the two contributions can be manipulated 
with the aid of the recursion relations between spherical 
Bessel functions: 


21 + 1 


ji{x) = ji-i{x) + ji+i{x), 


(21 + - (1 + ^)ji+i{x), 


to get, with the definitions JO(r) = f^J^(r) 

JL(r)=J^(r)-JC(r), 


(D14) 


f and 


<^10 

UJ 

^ 1—1 




27ry3 
1 


271^3 


/ASJLM + 

(D15) 


3. Cas© cit] 


As in Sec. ID 11 we split the two contributions: 


^-0 I 2 

^Im ^Im ' ^Im 


For 1 = 0, and recalling that Fqo = l/V^: 


J=o _ 
Clm — 


47r 


dp Wim(p) 


1 


(D19) 


1 


d'^r J,^(r)^=jo(fcc), 


1 

7W 


dp [pyim(p)]^ d^r 3a;{r)jo{kr). 


(D20) 

The result of the integrals in the shaded area above is: 


m = 1 : 2W- 


m = 0 : 2\ — 
3 


m = —1 ■. 2\ — 
3 

(D21) 


With which we reach: 


2. Case coo 

In the J = 0 case the contribution corresponding to 
I = j — 1 = —1 does not exist (see App. [B|), so the only 
contribution comes from 1 = 1 : 


|coo= ^ J dpWl^ip) J d^r 3^{r)Y{j^ji{kr) 

7n— — l 

rri—1 p p 

= y '^P (p)^ Fi^fi;(p) J d^r 3^{r)Y*^ji{kr). 


^10 

-^ 1-1 


For I = 2: 




J =2 _ 
^Im ~ 

—47r 


7n—2 


r -\/6 . 


d^r J,^(r)jo(fcc). (D22) 


dp [pYimip)]^ Y2rn J d^Y 3^{r)Y*:fj^j2{kr), 

(D23) 


(D16) the shaded momentum space integrals contain triple 
products of spherical harmonics and can be solved using 
The integrals in the shaded area are conveniently Eq. ( 1061 ) . They result in five vectors for each m case, 

solved using Eq. (jC5p and the orthonormality of the which we list here as row vectors. From top to bottom, 

spherical harmonics. After the sum in m we get: the row vectors corresponds to m = 2,1, 0, —I, —2: 
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TO = 1 —>■ 


TO = 0 —>■ 


TO = —1 —>■ 


-1 

7^ 


-1 

7^ 


-1 

7l5 


0 0 76 ) 
0 -73 0 ) 
10 0 ) 
0 0 0 ) 
0 0 0 ) 

0 0 0 ) 
0 0 73 ) 

0-20) 

73 0 0 ) 

0 0 0 ) 

0 0 0 ) 

0 0 0 ) 

0 0 1 ) 

0 -73 0 ) 

76 0 0 ) 


The following result is reached after taking steps par¬ 
allel to those taken in Sec. ID II for 
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d^Y js 7j^(r) f-J<^(r)| j 2 (/cr). 
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